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Abstract Optical and radio afterglows arising from the shocks by relativistic 
conical ejecta running into pre-burst massive stellar winds are revisited. Under 
the homogeneous thin-shell approximation and a realistic treatment for the lateral 
expansion of jets, our results show that a notable break exists in the optical light 
curve in most cases we calculated in which the physical parameters are varied 
within reasonable ranges. For a relatively tenuous wind which cannot decelerate 
the relativistic jet to cause a light curve break within days, the wind termination 
shock due to the ram pressure of the surrounding medium occurs at a small radius, 
namely, a few times 10^^ cm. In such a structured wind environment, the jet will 



X 



pass through the wind within several hours and run into the outer uniform dense 
medium. The resulting optical light curve flattens with a shallower drop after 
the jet encounters the uniform medium, and then declines deeply, triggered by 
runaway lateral expansion. 
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1 INTRODUCTION 

Gamma-ray Bursts (GRBs) are becoming gradually understood since the BeppoSAK af- 
terglow era (Piran 1999; van Paradijs, Kouveliotou & Wijers 2000; Cheng & Lu 2001; 
Meszaros 2002). Standard afterglow models have been set up within the hydrodynamical 
context of relativistic external shocks of spherical explosions running into either inter- 
stellar medium (Sari, Piran & Narayan 1998) or stellar winds (Chevalier & Li 2000), 
with synchrotron emission as the main radiation mechanism. Granot & Sari (2002) have 
obtained more accurate results by using the well-known Blandford & McKee (1976) self- 
similar inner structure of relativistic blastwave instead of the thin shell approximation. 
★ E-mail: xfwu@nju.edu.cn 
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However, the famous energy crisis under isotropic assumption (Kulkarni et al. 1999), 
together with the sharp decHne of some well observed optical afterglows led to the con- 
jecture of a jet-like form of this phenomenon (Rhoads 1997; Sari, Piran & Halpern 1999). 
Because of the relativistic beaming effect, the early afterglow from a jet is no different 
from a spherical blastwave. As the jet decelerates and the beaming effect weakens, the 
whole surface of the jet unfolds to the observer and the light curve declines more deeply 
with the decreasing radiating extent compared to the spherical blastwave. Additionally, 
the jet will decelerate much more rapidly when a significant or runaway lateral expansion 
takes place. Based on the homogeneous thin shell assumption, many authors have nu- 
merically studied jet dynamics and the behavior of the afterglow light curve (Panaitescu, 
Meszaros & Rees 1998; Huang et al. 2000a; Huang, Dai & Lu 2000b; Modcrski, Sikora & 
Bulik 2000). Tremendous efforts have been devoted to the fitting of the GRB afterglows 
of interest (Panaitescu & Kumar 2001, 2002). The jet kinetic energy is found to be sur- 
prisingly similar in different bursts, being tightly clustered around 3 x 10^" erg. Frail et 
al. (2001) found the genuine gamma-ray energy releases are also clustered around 5 x 10^° 
erg after the jet initial aperture 6j inferred from afterglows are accounted for (a factor of 

3 is amplified by Bloom, Frail & Kulkarni 2003). The total energy budget of a GRB is 
therefore close to that of a common supernova. Late time X-ray luminosity has been rec- 
ognized to be largely independent of the types of environment (Kumar 2000; Freedman 
& Waxman 2001). Subsequent statistics of the X-ray luminosities of several GRBs at 
10 hours since bursts led to the conclusion that Lx,iohr corrected by the jet aperture is 
again clustered around lO'*'* — 10^^ erg s~^ (Berger, Kulkarni & Frail 2003a). The above 
three independent estimations all identify the GRB as a standard energy reservoir and 
hence as a possible probe to the Universe. 

The discovery of optical transient of GRB 030329 (Price et al. 2003) and its associa- 
tion with SN 2003dh (Stanck et al. 2003) provides strong evidence for GRB-Supernova 
connection and removes any lingering doubt on the association of GRB 980425 with SN 
1998bw. These two associations together with several GRBs with late time re-brightening 
strongly imply near simultaneity of GRBs and SNe, with trigger time difference no more 
than 1 to 2 days (Wu et al 2003; Hjorth ct al. 2003). Affhough the detailed physical 
processes arc still uncertain, the central engine of at least long GRBs has thus been 
confirmed to be the core collapse of massive stars, called collapsars by Woosley (1993) 
as well as hypcrnovae for their energetics by Paczyiiski (1998). Previous studies favored 
the interstellar medium (ISM) as the environment of GRBs, albeit there were indications 
of wind environment for several GRBs within the spherical wind interaction model (Dai 
& Lu 1998b; Meszaros, Rees & Wijers 1998; Chevaher & Li 1999, 2000; Li & Chevaher 
1999, 2001, 2003; Dai & Wu 2003). The preference is partly due to the conclusion, based 
on the previous works (Kumar & Panaitescu 2000; Gou et al. 2001), that a jet in a stel- 
lar wind cannot produce a sharp break in the light curve. In fact, radiative loss of jet 
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energy was improperly neglected in these works, because the radiative phase lasts for 
hours or even one day due to the large wind density in the early afterglow. This fast 
cooling radiation will reduce the jet energy by about one order of magnitude. Realistic 
jets with energy losses in the stellar wind environments arc found to be consistent with 
several GRB afterglows (Panaitescu & Kumar 2002). Another origin of the bias against 
the wind environment arises from the relatively large or even impossibility of fit in 
some bursts. This result may be understood if the complicated mass loss evolutions of 
the progenitors at different stages and their interactions with larger scale environments 
are taken into account (Chevalier 2003). However, the region near the progenitor (within 
sub-pc) will not be affected and will be reflected in the intra-day afterglow. Even in the 
case of a supernova taking place 2 days earlier than the associated GRB, the supernova 
ejecta moving with one tenth of the light speed will reach no more than 10^^ cm, which 
is the typical deceleration radius marking the beginning of the afterglows. 

According to the above review of the recent progress of the GRBs and the wind 
environment, it seems worthwhile to revisit the evolution and radiation of homogeneous 
jets in a stellar wind environment. In this paper, we will study the dynamics and radiation 
in the jet plus wind model in Section 2. Numerical results with afterglow light curves 
clarifying the effects of different parameters are presented in Section 3. In Section 4 we 
give our discussion and conclusions. 



The dynamics of a realistic jet with radiative loss and lateral expansion running into 
a homogeneous interstellar medium is well described in Huang, Dai & Lu (2000b) and 
Huang et al. (2000a). Gou et al. (2001) studied in detail the evolution of a jet in a stellar 
wind environment, but they did not consider the radiative loss. We will follow these works 
and consider the radiative loss which is especially important in the wind environment 
case. We will also extend our calculation of the light curves into the radio band. Here we 
first give a brief description of our improved model. 

2.1 Dynamics 

The evolution of the radius R, swept-up mass m, half-opening angle 6 and Lorentz factor 
7 of the beamed GRB ejecta with initial baryon loading Mej and Lorentz factor 70 is 
described by (Huang et al. 2000a; Huang, Dai & Lu 2000b) 



2 DYNAMICS AND RADIATION 
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dm Mej + em + 2(1 — e)7m' 



(4) 



where t is the observer's time, /? = \/ 1^ ~ ^/l^ ^ the proton number density of the 
surrounding medium, nip is the mass of a proton, and e is the radiative efhcieney. For a 
stellar wind environment, the number density is given by (Chevaher & Li 1999, 2000) 

n = AR-"^, (5) 

where A = M /AwmpVw = 3 x 10"^^ A* cm~^, with M being the mass loss rate and 

A _ M 1 
* lO-SMoyr-i ^ lO^km s"! ' ^ ' 

is the wind parameter. The comoving sound speed Cg is (Huang et al. 2000), 

e,^=^(^-l)(^-l)__l_^e^ (7) 

where 7 « (47 + l)/(37) is the adiabatic index, which is appropriate for both relativistic 
and non-relativistic equation of state (Dai, Huang & Lu 1999). 

To estimate the radiative efficiency e, we assume the shock-accelerated electrons and 
the amplified magnetic field in the ejecta comoving frame carry a constant fraction, 
and ^B, of the total thermal energy. The magnetic energy density and the minimum 
Lorentz factor of the shock accelerated power law electrons in the comoving frame are 
then determined by 

^=^B4^-(7-l)nmpC^ (8) 

OTT 7 — 1 

7. = a7-l)^^^^f^ + l, (9) 

mc(p~ 1) 

m,c being the electron mass and p, the index in the electron power law energy distribution. 
The radiative efficiency of the ejecta is determined by a combination of the available 
fraction of energy to radiation contained in the electrons and the efficiency of the radiation 
mechanisms (Dai et al. 1999), 

^ = ^e-73^^, (10) 

where = R/{'yc) and t'^y^ = 67rmec/((7T-B'^7m) are the comoving-frame expansion 
time and synchrotron cooling time, respectively. 



2.2 Synchrotron radiation and self-absorption 

The distribution of electrons newly accelerated by the shock is assumed to be a power 
law function of the electron kinetic energy. Recently Huang & Cheng (2003) stressed 
that the distribution function should take the following form, dN^/d'-jc oc (7c — 1)"^ with 
7m < 7c < 7m, where 7m — 10^(-B'/1G)^^/^. This is especially important in the deep 
Newtonian stage. Radiation will cool down the electrons and thus change the shape of 
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the distribution. The cooling effect is significant for electrons with Lorentz factors above 
the critical value (Sari, Piran & Narayan 1998), 

7c = (11) 

In the comoving frame, the synchrotron radiation power at frequency v' from electrons 
of known distribution is given by (Rybicki & Lightman 1979) 

where is the electron charge, v'^ = Ssiw&j^qeB' /AnrriQC is the typical emission fre- 
quency of the 7e electron, and 

POO 

F{x) =x K5/3{k)dk, (13) 

J X 

with K^j^{k) the Bessel function. As emphasized by Wijers & Galama (1999), the random 
pitch angle between the velocity of the electron and the magnetic field will have some 
effects on the modelling of GRB afterglows. We choose the isotropic distribution of i? 
and have 

= 51^. (14) 
16meC ^ ' 

The characteristic frequencies corresponding to 7c, 7m and 7m electrons are denoted by 

f^, v'^ and v'-^^. 

In the wind environment the early radio afterglow flux density is rediiced significantly 
by synchrotron self absorption (SSA). This effect can be calculated by the analytical 
expressions derived by Wu et al. (2003). We give the more direct and convenient formula 
for the optical depth by SSA in different electron distributions as follows: 

1. For 1 < 7c < 7m, 

dNi^{ Co(7e - 1)-' (7c < 7e < 7m), (1^) 

^^'^ \ Cl(7e - 1)-^^+^) (7m < 7e < 7m), 

where 

= [(^ - ^-r) + , ^ ,, (1 - |^"~!lj ]-^A^eie, (16) 

7c - 1 7m - 1 mm - 1) (7M - 1)^' 

= Co(7m - lY-\ (17) 

and A'eie is the total electron number of a jet element. The self- absorption optical depth 
is 

7c-'(^)-^/' < <l 

-6,^^-3 



7c-'(-7)"' K <^'< O, (18) 

7m'(;;;^)-(''+')/' « <i^'< ^'mi 



7M 
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where Cga is (Wu et al. 2003) 

c«a = 10.4(p+2)/(p + 2/3), 

while 



m 



(19) 



(20) 



27r(l-cos^)i?2mp 

is the column density through which the synchrotron photons will experience 
synchrotron self absorption before emerging from the jet surface. We do not consider 
here the corresponding correction on the optical depth by the exact distribution of 
electrons, since the SSA coefficients are deduced under the assumption of soft photons 
and ultra-relativistic electrons, and when the bulk of the electrons are in the 
non-relativistic region the emitting source has already been optically thin to SSA. 



2. For 7m < 7c < 7m, 

d7e 



C2ilo~l)-P (7m < 7c < 7c), 

C3(7e-l)-(^+'^ (7c<7e<7M), 



(21) 



where (Huang & Cheng 2003) 

C2=C3/(7c-l), 

.(7m-l)i-P-(7c-l)i-^' , (7c-l)-^-(7M-l)-^_i 



The optical depth is 



(7c-l)(p-l) 



+ 



f ^-(f+4)fZl^-5/3 



: Csa Cl X < 



7m 

7c~ 



(p+4)^_^^-(p+5)/2 



7M 



3. For 7c > 7m, we have 



(22) 
(23) 



(24) 



dA^^ 

d7e 

where (Huang & Cheng 2003) 



C4(7e-1)-P (7m<7e<7M), 



Ca 



The optical depth is 



gc ^ Ncoi 



ele 



7n 



-(P+4)/ ^ N-5/3 
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'^m 
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The radiation power is assumed to be isotropic in the comoving frame. Let G be 
the angle between the velocity of an emitting element and the line of sight and define 
fj, = cos 6, the observed flux density at frequency v from this emitting element is 

<? 1 + ^ 1 l-exp(-T[7(l-/?/x)(l + z);/]) flnV^ m .^^A t09i\ 

'^ = ,Hl-P,rAnDli..) r[jil-P,)il + z).] ^ ^1 " + ^H, (28) 

where the luminosity distance is 

dz' 



D^{z)={i+z)^ r 

•no Jo 



(29) 



with Om = 0.3, ilA — 0.7 and Hq = 65 km/s/Mpc adopted in our calculations. The total 
observed flux density can be integrated over the equal arrival time surface determined 

by 

iobs = (1 + z)t ={l + z) [ ^—-^dR = const. (30) 
J PC 

3 NUMERICAL RESULTS 

We investigate the effects of various parameters on the afterglow light curve in our realistic 
jet+wind model. For convenience, we chose a set of 'standard' initial parameters as 
follows: Eiso = 2 X 10^^ erg, 6*0 = 0.1, = 0, 70 = 300, = 0.1, = 0.1, p = 2.2, and 

= 1.0. ^0 is the jet half-opening angle while 0^ is the observer's viewing angle with 
respect to the jet axis. The GRB in our calculations is assumed to be at ^; = 1, and the 
luminosity distance is about Dl = 7.1 Gpc. 

Figure 1 shows the optical and radio light curves with Eiso varying between 10^^ and 
10^* erg while the other parameters are fixed at their standard values. The early optical 
light curves resemble roughly those of spherical ones (for analytical jet+wind model, see 
Livio & Waxman 2000). The theoretical temporal evolution of the optical flux density 
changes from t~^^'^ to t~(3p-2)/4 ^ ^-i i5 Jqj. p _ 2.2) when crosses the optical band, 
at time (ChevaUer & Li 2000) 

= 0.1(i±^)^/^(|^)^/-^(||)^/^Ei/^33-K'/' days, (31) 

where = Z//4.36 x 10^^ Hz is scaled to the i?— band frequency and where wc have 
adopted the convention Q = lO'^Qx- The actual indices of the light curves in our cal- 
culations deviate slightly from the ideal spherical predictions. Wei & Lu (2002a, 2002b) 
suggested that some sharp breaks come from a spectral origin. Nevertheless, the radiative 
loss of kinetic energy and the effect of equal arrival time surface (EATS) cause several 
differences in the temporal evolution. For the extremely high isotropic energy case, viz. 
-Eiso, 53 = 10, there are two breaks in the light curve, as it evolves from t^"-^''' to t^^-^, 
and finally approaches the jet-like behavior t^^ *^. This temporal behavior is also seen in 
the 'standard' -Eiso, 53 = 2 case, except that the scaling law evolves from to t~^'^ 

and finally reaches t^^ *^. For other low Eiso cases, the light curves have only one break, 
while before the jet break the temporal index /3 is ~ 1.3 — 1.4 and after the jet break 
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(} K, 2.0 [Fi, cx t''^). The theoretical change of /3 in the spherical wind model takes place 
at (see also, Chevalier & Li 2000) 

ie = 3.8 X lQ\l±lf{^fErl,,Atvl days, (32) 

which is much later than the moment when significant lateral expansion is included. After 

that Uc will cross the R band and (3 changes from (3p — 2)/4 to (3p — l)/4 (i.e. from 

1.15 to 1.4 for p = 2.2). The difference of [5 between the jet model and the spherical 

model is ascribed to the EATS effect (see also Fig. 2 of Kumar & Panaitescu 2000). The 

radiative loss of energy also plays an important role for this diversity of the light curves. 

A crude estimate of the energy loss can be made by synchrotron radiation within the 

spherical fireball model. When the fireball is in the fast cooling stage, the observed power 

is P® = — ^e— ^ and the fireball energy evolves as 
di® ie 

^iso = ^iso,i(^)-«% (33) 

where E[so,i is the initial energy and ti = tdcc = 1-5 x 10~^( — ^ — )Eiso,D3^^^ day 
is the initial time of the afterglow, which begins at the deceleration radius. The fast 
cooling stage ends at to — 2( — - — )(q~]")((Jy)^* days (Chevalier & Li 2000). The kinetic 
energy decreases to i?iso.o ^ 0.25ii'iso,i (^o = 0.1) at Iq. The synchrotron power in the 

slow cooling stage is P® ~ ^2_(^)-(p-2)/2;^iso _ -j>j^g fireball energy is 

3-p to i® 

Eiso,{ « i;iso,oexp[-- -] « 0.07Siso,i, (34) 

(3-p)(p-2) 

in which we have put p = 2.2 and = 0.1. Early works neglected such tremendous loss of 
the jet energy and assumed relatively large isotropic energy and jet opening angle, leading 
to inconspicuous jet breaks (Gou et al. 2001). However, we show that notable breaks of 
optical light curves exist in all cases, as can be seen in Fig. 1. Although the temporal 
indices both before and after the break (see Fig. 1, t ^scveralxlO* seconds) vary with 
i?iso, the change of the temporal index in each case is nearly the same, A/3 « 0.65, and 
is well within one decade of time. Kumar & Panaitescu (2000) found that in a stellar 
wind f3 increases by ^ 0.4 over two decades of time, while in a uniform density medium f3 
increases by ^ 0.7 within about one decade of time. However, their conclusion is based on 
ignoring the enca'gy losses by radiation. Our present results thus make the jet+wind model 
competitive with the jet+ISM model. Radio afterglows are affected by the synchrotron 
self-absorption due to the dense wind medium at early times. A rapid rise of radio flux 
density proportional to t^ in Fig. 1 is consistent with theoretical expectation. As Eiso 
decreases, the transition to the optical thin regime becomes later than tm, and the radio 
light curve changes from type D to type E of Chevalier & Li (2000). At late times, the 
radio afterglows decay as t^^ because of lateral expansion. 

The effect of the wind parameter A* on the afterglow light curves is shown in Fig. 2. 
The early optical light curves (less than ~ 10"^ s) and late-time radio light curves (when 
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jet lateral expansion is significant) show little difference in the three cases with different 
A^,. The jet break in the optical band is indistinctive when A^, = 0.3, contrary to denser 
wind cases, in which an obvious optical break occurs around one day. The early radio 
afterglow is suppressed more strongly in denser winds. The radio light curves changes 
from type D to type E as A^ increases. 

For a relatively tenuous wind which cannot decelerate the rclativistic jet to cause 
a sharp break around one day, the wind termination shock due to the ram pressure 
balanced by the surrounding medium occurs at a small radius, i.e. several times 10^^ cm 
(Ramirez- Ruiz et al. 2001). Different mass loss histories of the progenitors or different 
large-scale environments surrounding the winds will lead to diverse environments for 
the GRB afterglow (see the review by Chevalier 2003). At least four kinds of GRB 
environment have been considered so far, stellar wind (Dai & Lu 1998b; Meszaros et al. 
1998; Chevalier & Li 1999), the ISM, dense medium (Dai & Lu 1999) and density jumps 
(Dai & Lu 2002). The complex, stratified medium resulting from interaction between 
the winds or between a wind and an outer dense medium interaction has been recently 
proposed to explain GRB 030226 by Dai & Wu (2003). This complicated and more 
realistic environment has the potential of unifying the diverse media, and, especially, 
of explaining the peculiar afterglow of GRB 030329 that included a large flux increase 
and several fluctuations. Fluctuations subtracted from the power-law light curve of GRB 
021004 indicate clouds and shells around the wind of a Wolf-Rayet star, which is the 
progenitor of an SN Ib/c and a coUapsar (Schaefer et al. 2003; Mirabal et al. 2003). In 
the complicated wind case, the jet will reach the wind termination shock radius at the 
observed time (Ramirez-Ruiz et al. 2001; Chevalier & Li 2000) 

U = lM^)Al_rErUnout hrs, (35) 

where riout is the number density of the outer uniform medium in units of cm^"^. After 
this time the jet enters the outer uniform medium. We calculate the light curve of such a 
jet with A* =0.1 and riout = 1 cm^''. The density ratio of the outer mcdimn to the wind 
at the termination shock radius is ^ 4, whic;li will not lead to a reverse shock propagating 
into the jet. As illustrated in Fig. 3, the resulting optical light curve flattens from an 
initial t^^-^ before 6 hours to t^^'^ after entering the uniform medium, and then declines 
steeply as t^'^-^ after ^ 9 days due to nmaway lateral expansion. This result provides the 
flrst piece of evidence that in the tenuous wind case an obvious sharp break can be caused 
by the medium outside the wind termination shock. The radio flux density will increase 
by a factor of a few since tt and thereafter follows the behavior of the jet+ISM model, 
which exhibits a jet break and shows a late time flattening when the shock becomes 
spherical and enters the deep Newtonian phase. 

The previous estimate of the jet half-opening angle 9q comes from direct fitting of 
afterglows and observed jet break time (Panaitescu & Kumar 2001, 2002; Frail et al. 
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2001; Bloom et al. 2003). Statistics shows a large dispersion of ranging from 2° to 
~ 40°. Here wc calculate the effect of on the light curves, this result is shown in Fig. 
4. A large 6*0 reduces the jet edge effects. For the case of = 0.2, the temporal index 
/? for the R band varies from 1.1 to 1.55 and then to 2.0 at ~ 100 days. However, jets 
with relatively smaller Oq will experience the jet break at earlier times. For the cases of 
9o = 0.075 and 0.05, (3 evolves from 1.3 — 1.4, to 2.0 at around 1 day. The corresponding 
A/3 is 0.7 and 0.6, respectively. The radio light curves show no difference at early times, 
due to the relativistic beaming effect and the same energy per solid angle. 

As wc know, the strong correlation between Eiso and 9q makes the GRBs a standard 
candle. A reliable treatment of the effect of i?iso or Oq to the jet light curves should 
include their intrinsic connection. We calculate the afterglows under the assumption of 
standard energy reservoir, Ej « Eiso-^ = 2 x 10^^ erg, and illustrate the results in Fig. 5. 
The characteristic feature of sharp breaks still remains for the narrow jets. For the case 
of 00 = 0.2, (3 for the R band evolves from 1.16 to 1.7, which is relative shallower than 
typical jet-break. The temporal index will further approach 1.86 tens of days later. For 
a standard candle with smaller angle of Oq = 0.05 (^o = 0.025), /3 evolves from 1.25 (1.3) 
to 1.85 (1.9) and then reaches 2.05 (2.06). Comparing with Fig. 4, the light curves in Fig. 
5 are more diverse at early times because the energy per solid angle is not a constant. 
Larger 9o will result in a smaller value of energy per solid angle and cause a lower level of 
flux density of both the optical and radio afterglows. However, the late time light curves 
are almost the same. 

We also investigate the effects of and on the light curves. Their most important 
effect is in the jet dynamics due to radiative loss, as illustrated in Eqs. (33) and (34). 
Figure 6 shows the effects of on the afterglows. In the case of = 0.2, (3 changes 
from 1.2 to 1.95. In the case of = 0.4, (3 changes from 1.5 to 2.2. Both cases show 
sharp breaks of Af3 w 0.7 at aroimd one day. The effect of is more complicated, since 
changing ^b significantly alters the type of the optical light curve. According to Eq. (32), 
tc is very sensitive to ^b- A small ^b will result in a type; B optical afterglow, in which 
f3 = (3p — l)/4 (Chevalier & Li 2000). The reason is that the three characteristic time 
scales depend on ^b as ^g, to ^b and oc ^g''''. Decreasing ^b will alter these 

three times significantly and lead to a type B afterglow, i.e. tc < f-a < ^m- Figure 7 shows 
that, for ^b = 10^^ and 10^'^, the f3 of the optical light curve evolves from 1.6 at early 
times to ^ 2.1 ten days later. Although A/3 « 0.5, the break is still obvious since it is 
completed within one decade of time in these cases. 

It is interesting to study the effect of the electron energy index, p, on the light curves. 
Theoretically p is very likely to lie within 2.2 - 2.4, but observations of GRB afterglows 
indicate that p may cover a rather wide range. Figure 8 shows the behavior of the 
afterglows with different p. It should be pointed out that the jet dynamics is not affected 
by p in our considerations. The light curves show no difference at very early times. For 
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the optical afterglow, /? evolves from 1.28 to 2.0 and 2.25 at very late times in the case of 
p = 2.4, while from 1.48 (1.6) to 2.2 (2.4) and finally to 2.55 (2.8) in the case oi p = 2.6 
{p = 2.8). The former breaks with A/3 ~ 0.7 - 0.8 are sharp enough to make the jet + 
wind model capable of fitting most of the observed breaks in afterglows. 

It is also important to determine the jet initial Lorentz factor through fitting the 
afterglows. Previous works have discussed the optical flash arising from a reverse shock 
when a fireball shell interacts with its circum-progenitor wind (Chevalier & Li 2000; Wu 
et al. 2003). The Lorentz factor can be determined from the optical flash. The imprint 
of Lorentz factor is also expected to appear in very early afterglows. Figure 9 illustrates 
the difference in the early optical and radio afterglows caused by different initial Lorentz 
factors. We see that lower Lorentz factors give lower flux densities at early times. However, 
late time afterglows depend mainly on the total energy of the jet and so the light curves 
in Figure 9 differ from each other only slightly at late stages. 

Lastly, we examine the effects of the viewing angle within our jet + wind model and 
illustrate the results in Fig. 10. This study may be of some help to our understanding of 
orphan afterglows, which are afterglows whose parent bursts are not observed because 
they lay off-axis with respect to our line of sight (Huang, Dai & Lu 2002). 

The actual profile of GRB jets is likely to be Gaussian (Zhang & Meszaros 2002a). 
It is interesting that the afterglow behavior of a Gaussian jet is very similar to that 
of a simple uniform jet, especially when the uniform jet is assumed to have no lateral 
expansion. This conclusion is supported by the hydrodynamical simulations of Kumar 
& Granot (2003), who showed that the lateral expansion velocity (ve) of a Gaussian jet 
is significantly less than the local sound speed. In other words, we can approximate a 
Gaussian jet with a uniform jet with lateral expansion velocity set at ve = 0. In Fig. 
11, we illustrate the effect of vg on the afterglow light curves. It is clearly shown that 
although the absolute intensities are different, the jet breaks in the i?-band light curves 
are similar in the both cases. In the vq = case, /? evolves from 1.0 to 1.7, so that the 
light curve breaks by about A/3 = 0.7 at around one day. 

4 DISCUSSION AND CONCLUSIONS 

The recently observed association of GRB 030329 and SN 2003dh strongly suggests a 
massive star origin for long GRBs. In the context of massive star evolution, a type Ib/c 
SN is expected to explode from a Wolf-Rayet star. The environment is thus likely to be 
a high speed stellar wind ejected from the progenitor. Despite the above basic physical 
reasoning, results of previous works favored instead an interstellar medium environment. 
This preference was partly due to the fact that the previous authors thought that the 
jet+wind model was unable to produce sharp breaks in the optical afterglow light curves. 
In this paper we have revisited a realistic jet+wind model and studied the effects of 
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various parameters, such as E^iso, ^o, ^e, Cb, etc.. Our more realistic model includes 
radiative energy loss and lateral expansion as well as the equal-arrival time surface effect. 
Inverse Compton scattering of synchrotron photons is not considered in this work, which 
can additionally decrease the cooling electron Lorentz factor 7c, leading to a further 
decrease of the jet energy. We find that obvious breaks can in fact be seen in the optical 
light curves. Temporal evolution of the jet energy due to radiative loss may be the main 
ingredient that gives rise to the notable break. Inverse Compton spectra components 
emerge always in the X-ray band, and will not significantly affect the temporal behavior 
of the optical afterglows in our calculations. Our results strongly suggest that breaks in 
the light curve can also be produced by jets expanding into stellar winds. The change of 
temporal index A/3 of these breaks, generally ranging from ~ 0.6 to 0.8, happens from 
less than a day to several days since the GRB while the transitions are generally well 
within one decade of time. The smallest break A/3 « 0.5 in our calculation occurs in a 
very weak wind. In reality, the stellar wind may be surrounded by an outer homogeneous 
medium, which would greatly complicate the physics of afterglows. 

Despite the overwhelming success of this simple jet model, the structure of GRB jets 
has recently been intrigued to the model constructors, although it has already been con- 
sidered in early works (Meszaros et al. 1998; Dai & Gou 2001). The motivation is driven 
by the large dispersion of both the isotropic gamma-ray energy releases and the jet aper- 
tures, compared to their tight correlation resulting in a standard candle as discussed 
above. There are two simple treatments of the jet structure that keep to the axial sym- 
metry. One common treatment is to assume the energy per solid angle and the Lorentz 
factor (therefore the baryon loading) decrease as power law functions of the angle from 
the jet axis. There exists a uniform and very narrow inner cone at the center in order to 
keep the total energy not finite. This kind of structured jet has the advantage of being 
able to explain the observed larger jet aperture with lower isotropic luminosity simply 
as the viewing angle effect, if the power law index of the energy distribution is about 
—2. It is also capable of reproducing the sharp breaks of some optical afterglows, as well 
as the observed luminosity function (Rossi, Lazzati & Rc!es 2002; Zhang & Meszaros 
2002a). As the transverse gradient of energy density is much smaller in the structured 
jet than in the homogeneous jet which has a definite boundary with the environment, 
the lateral expansion almost never approaches the local sound speed in the co-moving 
frame and can essentially be neglected in the analytic solutions (Kumar & Granot 2003). 
Detailed reexamination of structured jet is now available both numerically (Kumar & 
Granot 2003) and analytically (Wei & Jin 2003; Granot & Kumar 2003; Panaitescu & 
Kumar 2003; Salmonson 2003). However, the most probable profile of power law struc- 
tured jet is the one with the energy index —2 and a constant Lorentz factor, constrained 
by the existing afterglows (Granot & Kumar 2003). The other treatment is to assume 
the energy per solid angle as a Gaussian function of the angle from the jet axis (Zhang 
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& Meszaros 2002a; Zhang et al. 2004). This Gaussian jet, despite relatively small lat- 
eral expansion, will not deviate from the simple homogeneous jet significantly in their 
afterglow behaviors (Kumar & Granot 2003; Zhang & Meszaros 2002a). 

Although the power law structured jet model has the potential virtue of interpreting 
the large dispersion of GRBs within a unified picture, it is still confronted by several 
difficulties. (1) It is not supported by the simulations of rclativistic jet formation during 
the core collapsing of massive stars by MacFadyen & Woosley (2001). In their J32 model, 
the profile of the energy distribution at the emergence from the progenitor envelope is 
much better described by a Gaussian function, i.e. a Gaussian jet. In fact, the index 
should be —4 rather than —2 between 5° and 10° if a power law energy distribution is 
assumed. Zhang, Woosley & Heger (2004) have calculated the propagation of a rclativistic 
jet within its massive stellar progenitor (see also Zhang, Woosley & MacFadyen 2003). 
The ultimate structure at emergence is characterized by a uniform core component of a 
higher Lorentz factor spanning a few degree, with a sharp decline boundary adjoining a 
wider component with a lower Lorentz factor. (2) The power law structured jet model 
unifies most GRBs at the expense of introducing another adjustable parameter, 9c, the 
half-opening angle of the uniform core component. Different central engine properties 
will lead to different jet half-opening angle 9 (MacFadyen & Woosley 2001). To maintain 
the merit of power law structured jet model we have to regard 9c as a universal constant 
for all GRBs, and this has to be confirmed by observations through data fitting. (3) The 
most probable distribution of the Lorentz factor is a homogeneous one within the jet. 
But it is difficult to imagine such a distribution can be sustained at large angles and not 
affected by abundant baryon loading outside a very narrow funnel. The homogeneous 
jet model also appreciably suffers from this problem for a few GRB afterglows with 
large apertures, and can be settled by considering a two-component model as recently 
proposed (Berger et al 2003b; Huang, ct al. 2004). (4) Another difficulty for the power 
law structured jet model may be the prc-brcak flattening predicted by a few authors 
(Wei & Jin 2003; Kumar & Granot 2003; Granot & Kumar 2003). The flattening is due 
to the emergence of the core component before the jet break of tlic; light curve whilc^ the 
observer's line of sight is outside 9c. Such a flattening will even develop to a prc-brcak 
bump if 9c is within 1° (Salmonson 2003). This new type of bump has been proposed 
to explain the peculiar behavior of GRB 000301C, as an alternative to energy injection, 
density jump and microlensing scenarios (Dai & Lu 1998a; Zhang & Meszaros 2002b; 
Dai & Lu 2002; Garnavich, Loeb & Stanek 2000). However, it is suspect because some 
well-observed afterglows with large jet aperture such as GRB 000926 [9-^ = 8.1°) have not 
shown this pre-break bump behavior (Panaitescu & Kumar 2002). Future observations of 
very early afterglows in the upcoming Swift era would help to discriminate the structure 
of GRB jets. As proposed by Zhang & Meszaros (2002a), actual profiles of GRB jets may 
be Gaussian. The sideways expansion can be neglected since the gradients of physical 
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parameters in a Gaussian jet are much smaller than in an ideal uniform jet with a clear- 
cut lateral boundary to its environment. In such a case, a sharp break exists as well. With 
the self-absorption of synchrotron radiation included in our realistic jet+wind model, we 
are able to do broad-band fittings of observed GRB afterglows. 
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Fig. 1 Effects of the parameter Eiso on the optical (i?— band, upper panel) 
and radio (4.86 GHz, lower panel) light curves. The solid line corresponds to a 
'standard' jet with Eiso = 2x 10^^ erg, 0q = 0.1, 70 = 300, Cc = 0.1, = 0.1, 
and p = 2.2 running into a stellar wind with A* — 1.0. Other lines are drawn 
with only Eiso changed. 





Fig. 2 Effect of the parameter A* on the optical (i?— band, upper panel) and 
radio (4.86 GHz, lower panel) light curves. The solid line corresponds to a 'stan- 
dard' jet with Eiso = 2 X fO^^ ^^.g^ ^ q.I, 70 = 300, Co = O.f , £_b 0.1, and 
p — 2.2 running into a stellar wind with A* — 1.0. Other lines are drawn with 
only changed. 





Fig. 3 Optical (i?— band, upper panel) and radio (4.86 GHz, lower panel) after- 
glows from jets in a wind environment surrounded by a uniform medium. The 
solid line corresponds to a jet with Eiso ~ 2 x 10^^ erg, 9o — 0.1, jo — 300, 
= 0.1, = 0.1, and p = 2.2 running into a circum-stellar wind with 
A^^ = 0.1 and then entering the outer uniform medium of rtout = 1-0 cm^"^ 
when i? > i?c = 4 X lO" cm. An entire stellar wind situation is calculated 
(dotted line) for comparison. 





Fig. 4 Effects of the parameter 6*0 on the optical (i?— band, upper panel) and 
radio (4.86 GHz, lower panel) light curves. The solid line corresponds to a 'stan- 
dard' jet with Eiso = 2 X fO^^ ^^.g^ ^ q.I, 70 = 300, Co = O.f , £_b 0.1, and 
p — 2.2 running into a stellar wind with A* — 1.0. Other lines are drawn with 
only 00 changed. 





Fig. 5 Optical (i?— band, upper panel) and radio (4.86 GHz, lower panel) af- 
terglows from jets under the standard energy reservoir assumption, i.e. Ej = 

02 

Eiso-^ = 2 X 10^^ erg. The solid line corresponds to a 'standard' jet with 
Ei,o - 2 X 10^3 erg, = 0.1, 70 = 300, Co = 0-1, Cb = 0-1, and p = 2.2 
running into a stellar wind with = 1.0. Other lines are drawn with only 
changed while keeping Ej the same. 




Fig. 6 Effects of the parameter on the optical (_R— band, upper panel) and 
radio (4.86 GHz, lower panel) light curves. The solid line corresponds to a 'stan- 
dard' jet with Eiso = 2 X lO^^ ^^.g^ ^ q.I, 70 = 300, Co = 0.1, £_b 0.1, and 
p ~ 2.2 running into a stellar wind with — 1.0. Other lines are drawn with 
only Ce changed. 
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Fig. 7 Effects of the parameter on the optical (_R— band, upper panel) and 
radio (4.86 GHz, lower panel) light curves. The solid line corresponds to a 'stan- 
dard' jet with Eiso = 2 X 10^^ erg, 6*0 = 0.1, 70 = 300, Co = 0.1, £_b 0.1, and 
p ~ 2.2 running into a stellar wind with — 1.0. Other lines are drawn with 
only Cb changed. 




Fig. 8 Effects of the parameter p on the optical (_R— band, upper panel) and 
radio (4.86 GHz, lower panel) light curves. The solid line corresponds to a 'stan- 
dard' jet with E;so = 2 X f053 ^^.g^ ^ g.f , 70 = 300, Co = O.f , O.f , and 
p ~ 2.2 running into a stellar wind with — 1.0. Other lines are drawn with 
only p changed. 





Fig. 9 Effects of the parameter 70 on the optical (_R— band, upper panel) and 
radio (4.86 GHz, lower panel) light curves. The solid line corresponds to a 'stan- 
dard' jet with Eiso = 2 X 10^^ erg, 6*0 = 0.1, 70 = 300, Co = 0.1, £_b 0.1, and 
p ~ 2.2 running into a stellar wind with — 1.0. Other lines are drawn with 
only 7o changed. 




Fig. 10 Effects of the parameter 9^ on the optical (_R— band, upper panel) 
and radio (4.86 GHz, lower panel) light curves. The solid line corresponds to a 
'standard' jet with Eiso = 2 x fO^^ gj.g^ ^ q.I, 70 = 300, Cc = 0.1, = 0.1, 
and p = 2.2 running into a stellar wind with ~ 1.0. Other lines are drawn 
with only 9^ changed. 




Fig. 11 Effects of lateral expansion velocity vg on the optical (_R— band, upper 
panel) and radio (4.86 GHz, lower panel) light curves. The solid line corresponds 
to a 'standard' jet with E'iso = 2x10^^ gj.g^ ^ 0.1, 70 = 300, Co = 0.1, Cb = 0.1, 
and p — 2.2 running into a stellar wind with ~ 1.0, with vg — Cg. The dotted 
line corresponds to no lateral expansion, i.e. ve = 0. 



